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Nothing contained in any API publication is to be construed as granting any right, by 
implication or otherwise, for the manufacture, sale, or use in connection with any method, 
apparatus, or product covered by letters patent nor as indemnifying anyone from or against 
any liability for infringement of letters patent. 

This publication may be used by anyone desiring to do so. The Institute hereby expressly 
disclaims any liability or responsibility for loss or damage resulting from its use; for the 
violation of any federal, state, or municipal regulation with which an API publication 
may conflict; or for the infringement of any patent resulting from the use of an API 
publication. Every effort has been made by the Institute to assure the accuracy and 
reliability of the data presented. 
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FOREWORD 


This publication provides tables to correct hydrocarbon volumes metered under pressure 
to corresponding volumes at the equilibrium pressure for the metered temperature. The 
parallel publication in customary units is the Manual of Petroleum Measurement Stan- 
dards, Chapter 11.2.2. 

The table presented in this volume is also available from API as a computer tape, 
along with a manual containing the text information in this publication. 

Suggested revisions are invited and should be submitted to the director, Measurement 
Coordination Department, American Petroleum Institute, 1220 L Street, N.W., Wash- 
ington, D.C. 20005. 
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Chapter 11—Physical Properties Data 


SECTION 2—VOLUME CORRECTION FACTORS FOR METER PROVING AND 
HYDROCARBON COMPRESSIBILITY 


11.2.2M Compressibility Factors for 
Hydrocarbons: 350-637 
Kilograms per Cubic Metre 
Density (15°C) and — 46°C to 
60°C Metering Temperature 


11.2.2.1M SCOPE 


The purpose of this standard is to correct hydrocarbon 
volumes metered under pressure to the corresponding vol- 
umes at the equilibrium pressure for the metered tempera- 
ture. This standard contains compressibility factors related 
to the meter temperature and density at 15°C of the metered 
material. The corresponding customary version is Chapter 
11.2.2. 


11.2.2.2M HISTORY AND DEVELOPMENT 


The previous API standard for hydrocarbon compressi- 
bility, Standard 1101, Measurement of Petroleum Liquid 
Hydrocarbons by Positive Displacement Meter, was devel- 
oped from graphical correlations prepared in 1945. This 
standard was based on limited data with only a few points 
for pure fluids in the range from propane to pentane. No 
lighter mixtures and no effect of pressure on the compress- 
ibility factor were considered. In addition, no metric (SI 
version was available. 

In 1981, the Committee on Static Petroleum Measurement 
formed a subcommittee, the Hydrocarbon Compressibility 
-Group, to revise the compressibility tables of Standard 1101. 
As a result of an extensive literature survey, the data base 
found for the relative density portion of the table covers a 
broader range than that used in Standard 1101 but is lacking 
in data for unsaturated hydrocarbons. The data base was 
used to develop a mathematical model that includes the 
effect of pressure on the compressibility factor. The printed 
table produced from the model is the standard. This standard 
replaces the discontinued Standard 1101 and the first edition 
of Chapter 11.2.2, Compressibility Factors for Hydrocar- 
bons: 0.500-0.611 Relative Density Range and 20-128°F. 


11.2.2.3M TYPE OF STANDARD AND LIMITS 


The actual standard is the printed table of 251 pages that 
follows this text. The increments used in the table are 0.25°C 
and 2 kilograms per cubic metre. Interpolation to 1 kilogram 
per cubic metre in density is allowed. Compressibilities are 
in the usual units of reciprocal kilopascals but are calculated 
from two terms, A and B, and the pressure difference from 
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equilibrium, D,. This is necessary to obtain the desired 
accuracy in volume because of the important effect of pres- 
sure on the compressibility factor for light hydrocarbons. 
The range of the table is from —46°C to 60°C and from 
350 to 637 kilograms per cubic metre density (15°C), for 
use with pressure differences above equilibrium from 0 to 
15,200 kilopascals. 

The equation used to generate the table is given for those 
who wish to duplicate the table using their specific computer 
and language. Identical table information is available on a 
computer tape. The use of this computer tape to verify 
individually developed computer subroutines is highly rec- 
ommended. 


11.2.2.4M EXAMPLE USE OF THE STANDARD 


In this standard, the compressibility factor (F) is used in 
the normal manner for volume correction (* denotes mul- 
tiplication): 


Co = VelVmn = WML — F * Dy) 
Where: 


Cy = correction factor for pressure. 
V. = volume at the equilibrium (bubble point) 
pressure, P,. 
Vi, = volume at the meter pressure, P,,. 
DiS Pao Re 
P, and P,, may be in either kilopascals gage or kilopascals 
absolute, but both must be in the same units. 


As an example, calculate the volume at equilibrium pres- 
sure of 1000 cubic metres (V,,) of a material with a density 
(15°C) of 530.4 kilograms per cubic metre metered under 
a pressure of 5000 kilopascals at a temperature of 5.1°C. 
The equilibrium pressure (P,) for this material at 5.1°C is 
450 kilopascals. The rounded density and temperature val- 
ues of 530 kilograms per cubic metre and 5.0°C yield an A 
factor of 281,093 and a B factor of 5.504. The compress- 
ibility factor (F) is calculated as follows: 


F = 1A + D,*B) 
1/[281,093 + (5000 — 450) » 5.504] 
0000003267 


II 


The value for F is rounded to the ninth decimal place, to 
the maximum of four significant digits. ; 
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= L/[1.0 — 0.000003267 « (5000 — 450)] 
= 1.0151 
The value for C,, is rounded to the maximum of four decimal 
places. 
V. = Vin bod Co 
1000 « 1.0151 
= 1015.1 cubic metres 


The value for V, is rounded to the nearest 0.1 cubic metre. 


11.2.2.5M DATA BASE 


An initial 2278 data points were obtained from the lit- 
erature for pure fluid compounds and mixtures of light hy- 
drocarbon liquids. These data were examined to eliminate 
data for gases, data with large errors, and data with other 
abnormalities. The final data base used in this standard 
consists of 1724 data points from 13 sources (see Table 1). 
This metric standard was derived from the data base in U.S. 


customary units, so all discussion of the data base is limited @ 


to customary units. 

The ranges of the experimental data were relative den- 
sities (60°F/60°F) from 0.3477 to 0.6312, temperatures from 
— 28°F to 160°F, and pressure differences from 41 to 2036 
pounds per square inch gage (see Figure 1). The actual 
ranges for the standard, as determined by an API survey, 
are relative densities (60°F/60°F) from 0.350 to 0.637, tem- 
peratures from — 50°F to 140°F, and pressure differences 
from 0 to 2200 pounds per square inch gage. Hence, some 
portions of the standard represent extrapolated results. The 
uncertainty analysis presented in 11.2.2.7M may not be 
valid for these extrapolated portions. For the lower relative 
densities, 140°F is above the pseudocritical temperature at 
which liquid exists. For these fluids, the range is restricted 
to 96 percent of the pseudocritical temperature. 

The data set contains 46 different mixtures of normal 
hydrocarbons from methane to decane. The compositions 
of the mixtures are listed in Table 2. The use of the standard 
for compositions not close to those in the data base repre- 
sents an extrapolation whose results may have a greater 
uncertainty. 
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Table 1—Summary of Data Base 


Pressure 
(pounds 
per square 
inch page) 


Number 
of Data 
Points 


Relative 
Density Temperature 
(60°/60°F) (CF) 


32-140 
40-130 
70-160 
— 20-120 
— 20-100 
— 25-63 
— 28--66 
— 13-32 
32-140 
— 28-80 
— 28-100 
32-122 
122-149 


Sample References 


NGPA/TP2 
NGPA/TP1 
Cal Tech 


180-2000 

150-2000 21 
100-2000 9, 10, 13, 14, 15, 16, 17 
100-1500 1 

300-1600 

198-1788 

320-2200 

460-2100 

140-2120 

242—2040 

121-2130 

121-1477 

400-1465 


0.35-0.61 
0.35-0.51 
0.50-0.63 
Tulsa 0.35-0.51 
Manley/Swift 0.508 
Pope 0.356 
Straty 0.356 
Douslin 0,356 
Dittmar 0.508 
Haynes 0.51-0.58 
Ely 0.508 
Thomas 0.508 
-Teichmann 0.508 
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Figure 1—Limits of Data Base by Relative Density and Temperature 
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Table 2—Data Mixture Compositions (Mole Percent) 


iC; CG 


0.00 
5.74 
5.22 
4.78 
4.10 


3.59 
2.87 
2.26 
3.82 
3.47 


3.18 
2.73 
2.39 
1.91 
1.50 
1.97 
1.79 
1.62 
1.41 
1.23 
0.99 


OO AIN ARWNe 


0.77 0.37 
0.49 0.24 
0.00 0.00 
0.00 0.00 
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11.2.2.6M BASIC MODEL 


This metric standard is developed by converting metered temperature in degrees Celsius 
to degrees Rankine and density in kilograms per cubic metre at 15°C to relative density 
at 60°F. The converted parameters are then used in the compressibility factor equation 
from Chapter 11.2.2. In this metric standard, the compressibility factor (F) is used in 
the normal manner for volume correction: 


Cy = VolVm = UL = F « Dy) 


The effect of pressure on the compressibility factor (F') is not negligible for the range 
of conditions of the standard. Therefore, the pressure effect is included in the model by 
the use of two factors, A and B. The basic mathematical model used to develop this 
standard relates compressibility to relative density, temperature, and pressure difference 
as follows: 


F = 1A + Dy *B) 


Ax* 1075 = —2.1465891D—6 « TR? + 1.5774390D—5 « TR? « G? 
—1.0502139D —5 « TR? « G* + 2.8324481D—7 « TR? x G® 
—0.95495939D0 + 7,2900662D —8 « TR? « G? 
—2.7769343D —7 « TR? x G* + 0.03645838D0 «TR + G? 
—0.05110158D0 «TR *G + 0.00795529D0 «TR 
+9.13114910D0 «G 


A = A * 6.894757D0 


B «107-5 = —6,0357667D — 10 « TR? + 2,2112678D—6 « TR « G? 
+0.00088384D0 «G — 0.00204016D0 # G? 


D, = pressure above the equilibrium bubble point pressure, in consistent units of 
kilopascals gage or kilopascals absolute. 
DX = double-precision accuracy to the Xth power of 10. 
TR = temperature, in degrees Rankine. 
G = relative density (60°F/60°F). 


The use of higher powers of TR and G and other combinations of them did not improve 
the correlation. 
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11.2.2.7M UNCERTAINTY ANALYSIS 


The uncertainty in the compressibility factor is +10.8 
percent at the 95-percent confidence level. (The figure 10.8 
is 2.0 times the standard deviation of 5.4 percent, where 
2.0 is the two-tail probability value of a normal distribution 
for 1709 degrees of freedom at 95 percent.) These uncer- 
tainties represent the likelihood of the correlation’s ability 
to reproduce the data for a specific sample. They do not 
indicate how accurate the data are. In many cases, the ac- 
curacy of the experimental data is unknown. The corre- 
sponding uncertainty in volume and in C, is + 0.56 percent 
at the 95-percent confidence level, derived from a computed 
standard deviation of 0.28 percent, as described above. 

These volumetric uncertainties depend on operating con- 
ditions, the type of material, and the effect of pressure on 
the compressibility factor. They may not be true for the 
extrapolated portions of the standard. The regions where 
various uncertainties can be expected, averaged for all pres- 
sures, are plotted in Figure 2. The uncertainty for specific 
materials and temperature conditions can be obtained from 
this figure. For samples at lower relative densities, the un- 
certainty increases as the mixture’s critical temperature is 
approached. The correlation is valid for temperatures less 
than or equal to 96 percent of the pseudocritical temperature. 

If the effect of pressure on compressibility were ignored, 
there would be greater uncertainties in the volume. To il- 


lustrate this, the uncertainties in the calculated volume would 
range from 0.2 to 11 percent if a mean compressibility factor 
for 500 pounds per square inch, instead of the compressi- 
bility factor at the correct pressure, were used. This is from 
2 to more than 100 times the desired uncertainty of 0.1 
percent in the volume. Table 3 provides more details about 
the uncertainties due to ignoring the effect of pressure on 
compressibility at various conditions. 

In situations where either Chapter 11.2.2M or Chapter 
11.2.2 could be used to obtain corrected volumes, differ- 
ences in C,, can arise. Because of rounding, the increment 
in metered temperature of 0.25°C in Chapter 11.2.2M does 
not always yield C,, values equal to those from Chapter 
11.2.2, which uses 0.5°F increments. Table 4 shows the 
frequency of errors that can be expected when using tem- 
perature to the nearest 0.25°C, as opposed to the nearest 
0.5°F. In addition, maximum differences in C,, of +0.0001 
can be expected (at a frequency of 0.4 percent) because of 
conversion of pressure from pounds per square inch to kilo- 
pascals. It is therefore recommended that in cases where 
one party ordinarily uses metric units and the other party 
ordinarily uses customary units, the use of either Chapter 
11.2.2M or Chapter 11.2.2 should be agreed on before a 


transaction is made. 


Table 3—Effect of Pressure on Compressibility Factors 


Relative Temperature 
Density (°F) D, = 1000 psi 


0.360 -— 20 11.0 
68 26.0 


0.636 — 20 0.5 
140 5.1 


F Errors (%) 
D, = 2000 psi 


Cy Errors (%) 
D, = 1000 psi D, = 2000 psi 
27.0 1.5 : 
52.0 5.4 Hal 
1.5 0.4 0.2 
13.8 0.9 2:5 


Note: The values in the table are the errors when the compressibility factor at 500 psi, instead of the one at 


the measured pressure, is used. 


Table 4—Expected Frequency of Errors When Using 
Temperatures to the Nearest 0.25°C Versus the 
Nearest 0.5°F 


Absolute Difference in C,, 


0.0010 = Difference = 0.0015 
Difference = 0.0005 
Difference = 0.0001 


Frequency (%) 
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11.2.2.8M CALCULATION PROCEDURE 


This procedure is recommended for computers with 11 
or more significant digits. With some computers, this could 
require the use of double precision for all variables and 
constants, as shown in the following steps and example. 


11.2.2.8.1M Procedure Used to Generate 
Chapter 11.2.2M 


1. Vary temperature (TPC) from —46.00°C to 60.00°C 
(inclusive) in precise increments of 0.25; that is, 


TPC = —46.00, — 45.75, —45.50,..., 60.00 


2. Convert TPC to degrees Rankine (TR) and round TR 
to the nearest 0.01: 

TR = 1.80D0 « TPC + 32.00D0 + 459.70D0 

TR = INTCTR « 100.0D0 + 0.5D0) * 0.01D0 


The INT intrinsic function returns an integer by truncating 
all digits to the right of the decimal point. 

3. Vary density (DN) from 350 to 636 kilograms per cubic 
metre (inclusive) in precise increments of 2; that is, 


DN = 350, 352, 354, .. ., 636 


Extend the table to 637 kilograms per cubic metre. 
4. Convert the density in kilograms per cubic metre (DN) 
to density in grams per cubic centimetre (G): 


G = DN/1000.0D0 
5. Obtain three significant decimal figures by: 
G = INT(G * 1000,0D0 + 0.5D0) * 0.001D0 


6. Convert the density (15°C) in grams per cubic centi- 
metre to density (60°F) in grams per cubic centimetre: 

Gl =G 

G2 =Gs#G 

G3 = G2«G 

G4 = G2 «G2 


= 0.0D0 
G — 0.03689636D0 
= G + 1.24446244D0 » G1 


G — 0.63291568D0 « G2 
= G + 0.73861488D0 « G3 
G = G — 0.32478413D0 « G4 


7. Convert the calculated density (60°F) in grams per cubic 
centimetre to relative density using the density of water at 
60°F: 


Copyright por American Petroleum Institute 
Wed Feb 09 11:05:37 2005 


G = G/0.999012D0 
8. Obtain three significant decimal figures by: 
G = INT(G « 1000.0D0 + 0.5D0) « 0.001D0 


9. Calculate 96 percent of the critical temperature (TC): 


TC = 0.96D0 * (621.418D0 — 822.686D0 « G 
+ 1737.86D0 « G * G) 


10. If the temperature in degrees Rankine (TR) is more than 
the near critical limit (TC), void the calculation, output a 
message, and start the next density and temperature set. 
11. If the temperature in degrees Rankine (TR) is less than 
or equal to the near critical limit (TC), then calculate the 
A and B factors using Steps 12-18. 


12. Calculate the powers of temperature using temperature 
in degrees Rankine (TR): 

TR2 = TR«# TR 

TR3 = TR * TR2 


13, Calculate the powers of relative density: 


G2 =G«G 
G4 = G2 « G2 
G6 = G2 « G4 


14. Calculate the A factor: 
A = 0.0D0 
= A — 2.1465891D—6 * TR2 
= A + 1.5774390D—5 * TR2 « G2 
A — 1.0502139D—S * TR2 + G4 
+ 2.8324481D—7 « TR3 « G6 


0.95495939D0 
7.2900662D —8 « TR3 * G2 


2.7769343D —7 « TR3 « G4 
0.03645838D0 «TR » G2 
0.05110158D0 «TR «G 
0.00795529D0 «TR 
9.13114910D0 «G 


> Pr > > Pr >P> Pr Pr > 
oe en a a a ee 


15. Convert the A factor to kilopascals: 
A = A « 6.894757D0 


16. Round the A factor to the nearest whole number. 


A = INT(A « 100000.0D0 + 0.5D0) 


17. Calculate the B factor: 
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0.0D0 
B — 6.0357667D — 10 « TR2 


B + 2.2112678D—6 « TR « G2 
B + 0.00088384D0 «G 
B = B — 0.00204016D0 * G2 


www 
tl 


18. Round the B factor to the nearest 0.001: 
B = INT@ « 100000000.0D0 + 0.5D0) « 0.001D0 


19. Verify the table. Because of the complexity of the cal- 
culations, each term in the table should be verified against 
the standard. A tape of the table is available for use in 
verifying that the computer procedure will reproduce the 
standard. 


11.2.2.8.2M Use of the Procedure to Calculate 
Individual A and B Values 


Ll. Round the temperature to the nearest 0.25°C and ini- 
tialize the temperature in °C: 


TPC = XX.XX: —46.00 = TPC = 60.00 
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TPC is rounded to the nearest 0.25°C as follows: 
TT = INT(TPC): that is, truncation. 

DIFF = TPC — TT 

If DIFF = 0 then SIGN = 1.0 else SIGN = —1.0. 
DIFF = ABS(DIFF): that is, the absolute value. 


If DIFF < 0.125 then TPC = TT. 
If 0.125 = DIFF < 0.375 then TPC = TT 
+ 0.25 « SIGN. 


If 0.375 = DIFF < 0.625 then TPC = TT 
+ 0.50 * SIGN. 


If 0.625 = DIFF < 0.875 then TPC = TT 
+ 0,75 » SIGN. 


If DIFF = 0.875 then TPC = TT + 1.00 « SIGN. 


2. Round the density to the nearest 2 kilograms per cubic 
metre as follows: 


G = INT(G/2,0D0 + 0.5D0) «2 
3. Follow Steps 4—18 of the procedure given in 11.2.2.8.1M. 
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11.2.2.8.3M Example @ 


In this example, the temperature is 5.0°C and the density is 530 kilograms per cubic 
metre. 


NoTE: Results are shown to 18 digits for comparison purposes. Only 11 digits are required 
to generate the standard. Constants must use 1! digits, as indicated by the D exponent. 


1. Convert density at 15°C to relative density at 60°F: 


G = DN/1000.0D0 = 0.530000000000000000D + 00 
Gl =Ges 0.530000000000000000D + 00 
G2=GeG= 0.280900000000000000D + 00 


G2%*G = 0.148877000000000000D + 00 
G2 «G2 = 0.789048 1 00000000000D — 01 


0.0D0 = 0.000000000000000000D + 00 
G — 0.03689636D0 = — 0.368963600000000000D — 01 
G + 1.24446244D0 « GI = 0.622668733200000000D + 00 


G — 0.63291568D0 « G2 0.4448827 18688000000D + 00 
= G + 0.73861488D0 x G3 0.554845486177760000D + 00 
= G — 0.32478413D0 « G4 = 0.5292 18456 109094700D + 00 


G/0.999012D0 = 0.52974 1841048050173D + 00 
= INT(G « 1000.0D0 + 0.5D0) » 0.001D0 = —0.530000000000000000D + 00 


QAQAAAaAaAAAA 


2. Before rounding, the temperature is 5.00 = TPC. After rounding, the temperature 
is 0.500000000000000000D + 01. 
TC = 0.96D0 « (621.418D0 — 822.686D0 » G + 1737.86D0 * G * G) 


0.646616922240000000D + 03 
TR = 1.80D0 « TPC + 32.00D0 + 459.70D0 


0.500700000000000000D + 03 
3. Calculate the powers of temperature: 


TR2 = TR* TR = 0.250700490000000000D + 06 
TR3 = TR « TR2 = 0,125525735343000000D + 09 


4. Calculate the powers of relative density: 


G2 = GG = = 0.280900000000000000D + 00 
G4 = G2 « G2 = 0.789048100000000000D — 01 
G6 = G2 « G4 = 0.221643611290000000D — 01 


5. Calculate the A factor. 
Calculate the terms and sum to A: 
= 0.0D0 
A — 2.1465891D—6 « TR2 — 0.538150939 198659000D + 00 
= A + 1.5774390D—5 « TR2 « G2 0.572709488059854990D + 00 
— 1.0502139D-—5 « TR2 « G4 0.364961692917087106D + 00 
+ 2.8324481D—7 « TR3 « G6 0.115300476008567951D +01 


0.95495939D0 =  0.1980453700856795 10D +00 
+ 7.2900662D —8 « TR3 x G2 0.276853576564138604D + 01 


2.7769343D —7 « TR3 « G4 = —(0.180977793862079819D — 01 


MPMS LL.2.-2M-8b I | O?32e290 OOb9G1e ul 


SECTION 2—VoLUME CORRECTIONS FACTORS 11 


0.514584606164560798D + 01 
— 0.841503132453439202D + 01 
— 0.443181762153439202D + 01 

0.40769 1401465607976D + 00 

0.281093314409481085D + 01 


A + 0,03645838D0 «TR «G2 
= A — 0.05110158D0 «TR «G 
A + 0.00795529D0 « TR 

A + 9.13114910D0 *«G 

A = A « 6.894757D0 


> > > > 
I 


It 


Round and scale A to get the A factor rounded to the nearest whole number: 


A = INT(A « 100000.0D0 + 0.5D0) = 281093 


6. Calculate the B factor. 


Calculate the terms and sum to B: 


B = 0.0D0 

B = B — 6.0357667D—10 * TR2 = —0,151316966921568300D — 03 
B = B + 2.2112678D—6 *TR*G2 =  0.159690397175945700D — 03 
B = B + 0.00088384D0 «G =  0,628125597175945700D — 03 
B = B — 0.00204016D0 x G2 =  0,550446531759456999D — 04 


Round and scale B to get the B factor rounded to the nearest 0.001: 
B = INT(B  100000000.0D0 + 0.5D0) « 0.001D0 = 5.504 
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